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ABSTRACT 

Axial  flow  in  the  core  of  laminar  steady  trailing  vortex  from  the 
tip  of  a semi-infinite  lifting  wing  is  analyzed,  assuming  that  the  pressure 
gradient  is  determined  by  the  swirl  velocities  of  an  ideal  infinite  line 
vortex.  The  axial  and  lateral  variations  of  the  axial  velocity  depend  on 
the  strength  of  the  vortex  and  initial  axial  velocity  distribution  which 
must  be  specified  at  some  station  behind  the  wing. 
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I . INTRODUCTION 


Vortices  occur  naturally  and  are  also  of  technical  interest.  As  an 
example,  consider  the  trailing  vortex  which  is  formed  near  the  tip  of  a 


semi-infinite  lifting  wing  in  a uniform  steady  flow  u^.  The  vorticity 
shed  near  the  tip  of  the  wing  rapidly  rolls  up  into  a concentrated  core 
along  uq  and  the  flow  becomes  axisymmetric,  such  that  at  lateral  distance 
r -*•  00  the  velocity  u (r,z)  -*■  T /2irr  where  0 refers  to  the  angular 

u 0 

coordinate  and  z is  the  axial  distance  along  the  vortex  axis,  and  is  the 
circulation  around  the  vortex.  In  the  core  the  swirl  velocity  decreases  with 
downstream  distance  z.  The  pressure  is  constant  at  large  r and  hence  the 
pressure  increases  along  z and  in  the  core  there  must  he  an  axial  velocity 


difference  uz(r,z)  from  the  prevailing  velocity  uq.  These  simple  considerations 
show  the  existence  of  an  axial  velocity  difference  between  the  core  and  the 
surrounding  fluid.  The  theory  is  needed  to  find  its  detailed  distributions 


and  dependence  on  various  parameters.  This  paper  is  concerned  with  flows 
which  are  laminar  everywhere. 

Newman^  has  calculated  the  axial  velocity  variations  in  the  vortex  core 

by  assuming  that  both  u0  and  u^  are  small  compared  with  uq  and  the  equations 

2 

governing  u„  and  u become  uncoupled  and  may  be  solved  separately.  Dosanjh 

u 7, 

et  al  have  measured  the  velocities  in  the  turbulent  core  of  a vortex  and 


kali 


compared  the  results  with  Newman’s  calculation  for  the  laminar  case  assuming 
that  the  laminar  kinematic  voscosity  may  be  replaced  by  an  equivalent  tur- 
bulent exchange  coefficient. 

3 

Batchelor  has  analyzed  the  axial  velocity  in  the  vortex  core  which  is 

caused  by  the  axial  pressure  gradient  associated  with  the  decaying  vortex. 

In  effect,  he  assumes  that  u.  is  independent  of  u but  the  later  is  influenced 

0 z 


r 


by  the  former.  We  will  examine  at  the  end  of  this  paper  the  validity  of  this 
approximation.  Betchelor's  analysis  has  been  criticized  by  Tam  which  in  turn 
has  been  criticized  by  Herron^.  We  show  that  all  three  authors  have  missed 
some  essential  points.  We  give  here  the  correct  analysis. 

II.  GOVERNING  EQUATIONS 

The  equation  governing  uQ  in  trailing  vortex  under  boundary  layer 
(slender  vortex)  approximation  is  ^ 

3 , ^ % 2 3 , 2.  3 33  .V  (1) 

3z  o z 0 dr  r 0 3r  dr  r 

where  v is  the  kinematic  viscosity. 

All  other  investigators  have  invariably  assumed  that  the  axial 

velocity  difference  u (r,z)  between  the  core  and  the  surrounding,  and 

z 

the  radial  velocity  u^  (r,z)  may  be  neglected,  which  we  may  call  LINE 
VORTEX  approximation.  Presently  we  use  these  approximations  and  subsequently 
we  examine  the  limits  of  their  validity  and  estimate  the  effect  of  the 
neglected  terms  on  the  dynamics  of  the  vortex  flow.  For  the  line  vortex, 

Eq.  (1)  becomes 


u — u r 
o3z  0 


3_  3L  J 
V 3r  r 3r  r 


which  has  a self-similar  solution 


ufl(r,z)2irr  = T (l-e-n) 
o o 


where 


n = u r /4vz 
o 


Assuming  that  uf  is  given  by  Eq.  (3)  we  calculate  u in  the  following 
u z 

approximate  manner.  We  assume  that  the  radial  pressure  gradient  is  given 
approximately  by  the  equation. 


I = 

p dr  r 
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Using  Eqs.  (3)  - (5)  we  can  calculate  axial  pressure  gradient  , thus 

7 fz-  4>2&2  <p">' 


where 


p(n)  = /« 
n 


-t  2 

(1-e  V dt 


The  linearized  equation  of  motion  for  u is 


a i an  .a2  . i a . 

U r—  U T^+V  ( _ + — — ) U 

o 8z  z p 3z  . 2 r 3r  z 

3r 


III.  AXIAL  VELOCITY 


The  axial  velocity  difference  uz  la  determined  by  the  initial 

distribution  of  u and  by  the  pressure  gradient  along  the  vortex  axis 
z 

which  is  due  to  decaying  u in  the  vortex  core.  The  diffusion  Eq. 

e 

(8)  has  a unique  solution  for  the  sources  given  by  Eq.  (6)  and  an  inltla1 
condition  u*(r,ZQ)  prescribed  at  a station  zq.  Using  standard  Green 
Function  technique,  see  for  example  Luikov^,  the  solution  due  to  the  sources 
is  given  by 
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u (r,z) 
z 


(P(n)n) * e_ 


~2  x .2 
5 +8 


4(t-f) 


<t-t’) 
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88 


o'2(t-t ') 
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o 
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where  s = ru  /v  , t = u z/v  , ri  = s /4l  and  t = u z / \>  and  1 Is  the 
o o o o o o 

modified  Bessel  function. 

The  contribution  u*  (r,z)  due  to  the  initial  u*(r,zo>  is  again  obtained 
using  standard  methods  by  the  relation 


(9) 


(10) 


(11) 

In  the  present  notation  total  axial  velocity  is  the  sum  of  uz(r,z) 

and  u*(r,z)  . The  sources  given  by  (6)  have  a non-integrable  singularity 
z 

at  z = 0,  so  that  zq  cannot  be  zero.  The  initial  velocity  cannot  be  pre- 
scribed at  z * 0,  and  therefore  the  solution  is  not  independent  of  zq  . 

Unfortunately,  Eq.  (9)  cannot  be  integrated  in  terms  of  known  functions. 

We  have  numerically  integrated  it  for  z^u^/v  * 1 and  z/zq  = 10*  100,  and 
1,000.  We  have  also  integrated  it  for  another  Initial  station  zquq/v  = 10 

and  for  z/z  * 10,  and  100.  The  velocity  profile  u (r,z)/u  (o,z)  are 
o z z 

shown  in  Fig.  1.  The  maximum  non-demens ional  velocity  defect  uz(o,z)/uq, 
as  a function  of  z,  is  shown  in  Fig.  2 and  3. 
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for  z>>z 
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IV.  AXIAL  VELOCITY  USING  EQUIVALENT  SOURCES 


The  function  jp(n)rjj  of  the  source  distribution  is  shown  in  Fig.  4 
together  with  e 11  . The  area  under  the  two  curves  (i.e.,  the  integral  of 
sources  in  the  radial  plane)  is  the  same  and  the  curves  are  also  similar. 
Therefore,  we  may  approximate  Eq.  (8)  by 


- ~ (?£■) 
p V9z' 


r 0u 

. (-°)  2_£L_ 
*2ir  _ 2 

8vz 
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-r  u 
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4vz 


(12) 


Using  (12)  as  the  source  Eq.  (8)  has  the  solution 


2 

-r  u 


u (r,z)  _i  rn  ? 1 4vz 

— (— ) a — in  z/z  e 

u u 2ir  8vz  o 

o o 


(13) 


We  now  find  that 


2 

-r  u 


4vz 


u (r,z)/u  (o,z)  * e 
z z 


and 


(14) 


u (o,z)/u 
z o 


A)2 


1 

8vz 


in  z/z 

o 


(15) 


The  results  given  by  Eqs.  (14)  and  (15)  compare  favorably  with  the 
numerical  computation  as  shown  in  Fig.  1,  2 and  3,  respectively. 

We  note  that  u is  the  axial  velocity  defect  due  to  the  pressure 
gradient  in  the  vortex  core  starting  with  uniform  axial  velocity  at 
zq.  We  may  add  to  it  the  contribution  of  any  initial  axial  velocity 
excess  or  defect  using  Eq.  (10)  or  Eq.  (11)  as  appropriate. 


A brief  version  of  this  analysis  was  given  earlier. 


V.  DISCUSSION 


A solution  of  Eq.  (8)  for  the  sources  given  by  Eq.  (6)  was  given  by 

3 4 

Batchelor  which  was  criticized  for  non-uniqueness  by  Tam  which  in 

turn  was  criticized  by  Herron^.  The  difficulties  are  due  to  the  singularity 

of  the  source  distribution  at  z = 0,  see  Eq.  (6),  which  no  one  has  noticed. 

If  the  source  distribution  were  integrable  at  z = 0 then  it  might  have  been 

possible  to  get  a solution  which  does  not  depend  on  the  detailed  distribution 

of  sources  and  which  would  be  valid  for  large  z. 

We  can  estimate  the  neglect  of  radial  velocity  and  the  variations  in 

axial  velocity  in  Eq.  (2)  using  the  procedure  given  in  reference  6.  Integrating 

Eq.  (1)  we  have 


oo  !)  ? oo 

u / v”  u„r  dr  + / — uftu  r2dr  + 

o o 3z  0 o 3z  0 z 


"*or  = -vr  /it 


or  using  the  continuity  equation  and  F = 2irrun  , we  have 


u /"  r rdr  - /"  |-  (r  -r)u  rdr  - 2vr 

O O dz  0 0Z  0 Z 0 


convection 


diffusion 


The  ratio  of  the  neglected  terms  which  represent  the  radial  and  the  associ- 


ated axial  convection  of  angular  momentum  to  the  diffusion  term,  which  is 

retained,  can  be  calculated  from  the  analytical  solution  for  u based  on 

z 

equivalent  sources. 

Thus  we  find 

«b  a r 

/_  JT  (r  -r)u  rdr/-2vr  « (-J)  (r  /64irz  U ) (1* 

odzoz  ovo  oo 
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We  note  that  for  large  z the  axial  velocity  due  to  any  Initial  distribution 
does  not  contribute  to  the  above  ratio.  The  approximations  used  to  replace  the 
trailing  vortex  by  a line  vortex  are  valid  if  the  above  ratio  Is  much  less 
than  unity. 

The  comparison  of  this  theory  with  available  experiment  can  only 

be  made  when  Eq.  (18)  is  satisfied  and  the  flow  is  laminar.  We  do  not  know 

of  any  such  experiments.  In  general,  the  results  of  analysis  based  on 

laminar  flow  cannot  be  applied  to  the  flow  in  the  turbulent  state.  In  the 

case  at  hand,  it  is  not  possible  to  replace  the  kinematic  viscosity  by  a 

turbulent  exchange  coefficient  and  use  the  theoretical  results  presented 

here  for  comparison  with  experimental  results  for  turbulent  trailing 

vortices.  The  reasons  are  not  directly  relevant  here,  they  are  quite 

9 

extensive  and  are  presented  elsewhere. 
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Fig.  1.  Axial  velocity  defect  distribution  as  a function  n, 
e~n;  x z/zq=  1,000,  zq=  1;  A z/ Zq=  100,  zq  = 1 ; 

0 z/Zl  = 10,  zQ  = 1;  l"D  z/zq=  100,  zq  = 10,  + z/zq  = 10 

Fig.  2.  Axial  velocity  defect  distribution  on  the  vortex  axis  for 

z/z  from  1 to  7;  - approximate  solution  using  equivalent  sources; 
o 

x numerical  solution  using  exact  distribution  of  sources; 

Fig  3.  Axial  Velocity  defect  distribution  on  the  vortex  axis  for 

z/z  from  1 to  100;  - approximate  solution  using  equivalent  sources 
o 

x numerical  solution  using  exact  distribution  of  sources; 

Fig.  4.  Comparison  of  sources  with  equivalent  sources; 
x sources  distribuiton;  - equivalent  source  distribution. 
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